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The  most  promising  materials  that  can  be  used  as  anodes  in  next  generation  rechargeable  Li  batteries  are 
Sn  and  Si.  Upon  lithiation,  however,  both  Sn  and  Si  experience  a  300%  volume  expansion,  which  results 
in  significant  fracture,  and  therefore  their  commercial  use  is  inhibited.  Extensive  experimental  research 
has  yielded  that  embedding  or  attaching  Si  or  Sn  nanoparticles  in  a  carbon/graphite  matrix  diminishes 
their  mechanical  damage  and  allows  for  electrochemical  stability.  The  present  study  will  show  that 
linear  elasticity  can  predict  the  capacity  retention  of  such  nanocomposites  by  predicting  their  mechanical 
stability  upon  Li-insertion.  In  particular  (i)  a  previously  developed  theoretical  model  will  be  related 
to  experimental  observations  on  Si/sol-gel-graphite  nanocomposite  anodes,  (ii)  electron  microscopy 
images  will  be  presented  on  the  fracture  of  cycled  Sn02/C  nanopowders,  and  a  theoretical  model  will  be 
applied  to  predict  the  Sn02  particle  dimensions  that  will  limit  such  fracture. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Secondary  Li  batteries  are  the  most  rapidly  growing  high  energy 
density  power  sources;  their  applications  ranging  from  hybrid  and 
electric  vehicles,  to  biomedical  implantable  devices  [1].  Commer¬ 
cial  anodes  for  Li  batteries  comprise  of  graphite  since  it  offers  both 
electrochemical  and  mechanical  stability.  One  of  the  drawbacks, 
however,  is  that  graphite  has  a  low  Li  intercalation  (LixC6,  where 
0<x<  1),  and  therefore  provides  a  capacity  of  372  mAhg-1  [2]  upon 
lithiation. 

About  the  same  time  that  graphite  was  proposed  as  a  stable 
anode,  experimental  research  [3]  yielded  that  the  most  promis¬ 
ing  anodic  materials  are  Si  and  Sn,  which  allow  for  capacities  of 
4200  mAh  g-1  and  990  mAh  g-1,  respectively,  upon  the  formation 
of  Li  intermetallics  [4].  These  high  capacities  are  partly  attributed 
to  the  higher  Li-intercalation  in  these  materials,  namely  LixSn  and 
LixSi,  where  0  <x  <  4.4.  Despite  their  unique  electrochemical  prop¬ 
erties,  however,  Sn  and  Si  cannot  be  used  commercially,  as  their 
mechanical  behaviour  pays  the  penalty  of  their  high  reactivity  with 
Li.  In  particular,  during  maximum  Li  insertion,  Si  and  Sn  experience 
volume  expansions  of  300%  [4].  As  a  result  during  electrochemical 
cycling  (continuous  insertion  and  de-insertion  of  Li)  the  anodes 
expand  (lithiation)  and  contract  (de-lithiation)  leading  to  fracture. 
This  fracture  produces  active  material  that  is  no  longer  in  electri¬ 
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cal  contact  with  the  remainder  of  the  electrode,  which  therefore 
becomes  unable  to  respond  to  the  applied  voltages  necessary  to 
recharge  or  control  the  discharge  of  the  battery. 

The  capacity  of  bulk  Sn  and  Si,  therefore,  becomes  negligible 
after  the  first  few  cycles.  At  the  nanoscale,  however,  deforma¬ 
tion  mechanisms  are  less  severe,  and  better  capacity  retentions 
have  been  achieved.  In  [5]  it  was  illustrated  that  during  the  first 
Li-insertion  the  capacity  of  Sn  nanoparticles  (50-200  nm  in  diame¬ 
ter)  was  990  Ah  g-1 ,  but  it  dropped  to  21 0  mAh  g-1  after  20  cycles. 
Similar  trends  were  shown  for  Si  nanoclusters  that  gave  a  start¬ 
ing  capacity  of  2400  mAh  g-1 ,  which  reduced  to  525  mAh  g-1  after 
fifty  cycles  [6].  Post  mortem  transmission  electron  microscopy 
images  document  that  this  large  capacity  reduction  during  cycling 
is  attributed  to  the  severe  fracture  that  results  from  the  volume 
expansion  of  the  Sn  and  Si  nanopaticles. 

The  most  effective  way  to  constrain/accommodate  these  expan¬ 
sions  is  to  embed/attach  the  active  materials  with  respect  to  Li, 
in  a  less  active  matrix.1  The  matrix  constrains/buffers  the  volume 
expansions,  and  also  protects  the  active  site  surface  from  elec¬ 
trochemical  damage.  One  of  the  most  effective  such  matrices  is 
carbon/graphite,  since  once  maximum  Li-insertion  takes  place  in 
the  active  site,  additional  Li  can  be  intercalated  by  the  carbon. 

Niu  and  Lee  [7]  were  able  to  achieve  capacities  greater  than 
2400  mAh  g-1  for  20  cycles  by  embedding  Si  nanoparticles  in  a 


1  By  embedding  the  active  material  is  fully  surrounded  by  the  matrix,  whereas  by 
attaching  the  active  material  is  attached  on  the  surface  by  the  matrix  forming  an 
island  structure. 
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sol-gel-graphite  matrix,  while  Derrien  et  al.  [8]  were  able  to  attain  a 
capacity  of  500  mAh  g-1  for  1 00  cycles  by  embedding  Sn  nanoparti¬ 
cles  in  graphite.  In  addition,  to  embedding  the  nanoparticles,  it  was 
shown  that  attaching  Sn  particles  on  an  amorphous  carbon  surface, 
increased  the  carbon  capacity  by  85%  [9].  Instead  of  using  carbon  or 
graphite  as  the  matrix,  cellulose  fibers  have  also  been  successfully 
employed;  Si  nanoparticles  attached  on  cellulose  gave  capacities 
of  1300  mAh  g-1  for  30  cycles  [10],  and  Sn  attached  on  cellulose 
gave  a  capacity  of  500  mAh  g-1  for  40  cycles  [11].  Such  Sn  and 
Si  based  anode  materials  are  usually  referred  to  as  active/inactive 
composites,  since  the  Sn  and  Si  are  significantly  more  active  with 
respect  to  Li  as  compared  to  the  matrix.  It  is  clearly  seen  that  using 
such  nanocomposite  systems  as  anodes  results  in  a  much  better 
capacity  retention  than  using  Sn  or  Si  in  their  pure  form  [5,6].  A 
more  detailed  overview  on  various  configurations  and  materials 
chemistries  that  are  being  considered  as  next  generation  anodes  in 
Li  batteries  can  be  found  in  [1  ]. 

It  should  be  noted,  however,  that  not  only  materials  chem¬ 
istry,  but  also  the  particle  size  and  the  volume  fraction  of  the 
active  sites  affect  the  capacity.  In  [7]  for  example,  three  different 
Si/sol-gel-graphite  anodes  were  tested,  with  a  different  Si  content. 
The  capacity  retention  was  the  best  for  the  anode  with  the  lowest  Si 
content.  Furthermore,  the  promising  Sn/graphite  anode  presented 
in  [8]  had  a  mean  Sn  particle  diameter  of  35  nm,  whereas,  earlier 
studies,  that  used  similar  Sn/graphite  nanocomposites,  but  with  a 
Sn  average  diameter  greater  than  lOOnm  gave  very  poor  capacity 
retentions  [12]. 

Although  the  existing  experimental  evidence  can  provide  some 
insight  towards  selecting  appropriate  particles  size  and  volume 
fractions,  design  criteria  can  be  obtained  only  through  the  develop¬ 
ment  of  theoretical  models  that  can  predict  electrochemical  and/or 
mechanical  stability.  In  the  present  study,  therefore,  the  high  capac¬ 
ity  retention  that  can  be  obtained  for  active/inactive  anodes  will  be 
related  to  their  enhanced  mechanical  stability.  The  theoretical  pre¬ 
dictions  that  will  be  deduced  through  linear  elasticity  and  energy 
considerations  will  be  directly  used  to  interpret  existing  and  new 
experimental  data. 


2.  Mechanical  models  developed  for  Li  anode  fracture 

Although  there  exist  myriad  studies  on  the  experimental  aspects 
of  Li-electrodes,  there  is  limited  work  that  focuses  on  the  mechani¬ 
cal  issues  of  such  systems,  despite  the  experimental  indications  that 
these  systems  fail  due  to  severe  damage  and  fracture  [4,13].  Before, 
therefore,  Sn  and  Si  based  anodes  can  be  used  commercially  their 
fracture  mechanisms  must  be  understood. 

Wolfenstine  [14]  employed  linear  elasticity  to  predict  that  the 
particle  size  of  Sn  that  will  result  in  no  fracture,  during  lithiation,  is 
below  its  unit  cell  size.  Huggins  and  Nix  [15],  however,  predicted 
that  the  particle  size  that  will  inhibit  fracture  is  below  0.35  p,m  for 
brittle  materials  and  8  fxm  for  ductile  materials,  which  however,  is 
not  consistent  with  experimental  evidence.  This  discrepancy  from 
experiments  is  most  likely  due  to  the  fact  that  the  work  in  [15]  was 
based  on  the  fracture  analysis  model  of  a  thin  film  configuration 
and,  hence,  the  internal  stress  development  and  fracture  of  an  epi¬ 
taxial  thin  film  on  a  compliant  substrate  [16]  was  analyzed,  instead 
of  the  actual  radial  symmetry  that  is  present  in  active/inactive 
nanocomposites. 

The  mechanics  challenges  that  arise  for  electrodes  comprised 
of  active  nanoparticles  embedded  in  a  matrix  were  initially  under¬ 
taken  in  [17-19],  through  the  development  of  simple  physically 
based  linear  elasticity  and  fracture  mechanics  models.  The  goal  was 
to  establish  preliminary  easy-to-use  material/component  design 
criteria  for  nanostructured  anodes  [20]  with  optimum  capacity, 
strength,  and  electrochemical  stability.  In  particular,  under  con- 
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Fig.  1.  Experimental  data  exhibiting  the  capacity  fade  of  nano-Si  embedded  in  a 
sol-gel  graphite  matrix  (taken  from  [7]);  lower  Si  content  has  a  lower  capacity  fade. 


tinuous  electrochemical  cycling  (nano)  cracking  takes  place  at  the 
particle-matrix  interface  and  a  (nano)  damage  zone  is  developed 
as  documented  by  transmission  electron  microscopy  observations 
[  1 7,2 1  ].  On  assuming  that  this  damage  zone  is  comprised  by  a  num¬ 
ber  of  radial  cracks,  Aifantis  et  al.  [20]  were  able  to  derive  stability 
criteria  for  the  crack  growth  depending  on  the  geometry  assumed 
for  the  active  nanoparticles  (spheres,  cylinders  or  disks),  and  the 
associated  boundary  conditions. 

Strain  considerations  for  radial  geometries  for  the  active  sites 
have  also  been  recently  examined  in  [22],  which  looked  at  mod¬ 
erately  expanding  cathodes,  with  geometries  similar  to  those 
advanced  previously  in  [17-20],  but  also  accounted  for  stress- 
assisted  diffusion  similarly  as  in  [23].  A  more  recent  study  [24] 
highlights  how  stresses  are  built  up  in  porous  electrodes  which  are 
an  agglomerate  of  active  particles,  the  electrolyte,  and  the  inactive 
material. 


3.  Volume  fraction  considerations  in  Si-sol-gel-graphite 
(Si-SGG)  anodes 

In  order  to  understand  the  effect  that  the  Si  volume  fraction  has 
in  the  overall  capacity  retention,  Niu  and  Lee  [7]  considered  three 
different  anodes:  19.2  wt%Si-SGG,  9.8wt%Si-SGG  and  4.8wt%Si- 
SGG.  As  seen  in  Fig.  1  the  nanocomposite  with  the  highest  Si  content 
(19.2  wt%Si)  gave  a  higher  initial  capacity,  since  it  contained  more 
active  material,  but  after  25  cycles  the  capacity  decreased  by  18%. 
The  lowest  Si  content  (4.8  wt%Si)  anode  gave  a  lower  initial  capac¬ 
ity  which,  however,  remained  stable  for  25  cycles.  In  this  section  it 
will  be  shown  that  this  correlation  between  Si  content  and  capacity 
retention  is  attributed  to  the  underlying  mechanical  stability. 

In  [18]  linear  elasticity  and  axial  symmetry  was  employed  to 
study  the  fracture  of  nanocomposite  anodes  in  which  spherical 
active  sites  were  embedded  in  a  matrix;  a  unit  cell  is  shown 
in  Fig.  2.  The  fracture  that  occurs  during  electrochemical  cycling 
was  accounted  for  by  considering  a  damage  zone  at  the  active 
site/matrix  interface.  This  zone  was  characterized  by  multiple 
radial  cracks  of  length  p-a ,  whose  opening  tensile  stress  was  ctqq. 
Earlier  work  by  Dempsey  and  co-workers  [25,26]  concluded  that 
the  energy  release  rate  (G)  during  cracking  in  such  geometries  is 
given  by 

c  _  (1  -  VmYfpJ  r  (fr  +  2pXb  -  p)l 

2  nEmp3  b2  +  bp  +  p2 

where  a  and  b  are  the  radii  of  the  active  site  and  matrix,  p  is  the 
crack  radius,  vm  and  Em  are  the  Poisson’s  ratio  and  elastic  mod¬ 
ulus  of  the  matrix,  and  n  is  the  number  of  radial  cracks,  p  is  the 
internal  pressure  that  the  active  site  exerts  onto  the  matrix  upon 
its  300%  volume  expansion  during  maximum  Li-insertion  and  was 
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Fig.  2.  Configuration  of  unit  cell  used  in  analysis;  a  and  b  are  the  radii  of  the  active 
site  and  matrix,  A  is  the  free  expansion  of  the  active  site  if  it  were  not  constrained 
(by  the  matrix),  and  p  is  the  crack  radius. 


p  in  nm 


Fig.  3.  Griffith’s  criterion  for  various  volume  fractions  ( b  is  kept  1000  |xm),  of  Si 
nanospheres  (with  radius  a)  embedded  in  a  Y2O3  matrix.  Arrows  indicate  the  crack 
radius  at  which  fracture  will  stop. 

computed  as  [19] 

P  (1  —Pm)  3(1  —  vm)p2  (a  +  A)ES  \ 
a  2  2{b2  +  bp  +  p2)\  1-2 vs  J  ’ 

(2) 


where  vs  and  Es  are  the  Poisson’s  ratio  and  elastic  modulus  of  the 
active  site  and  A  corresponds  to  the  distance  the  active  site  would 
expand  to  if  it  was  not  constrained  by  the  matrix.  By  inserting, 
therefore,  the  material  parameters  for  Si  and  C  in  Eqs.  (1)  and  (2) 
the  energy  released  during  fracture  in  Si-C  nanocomposites  can  be 
predicted.  In  order  to  capture  the  effects  of  the  volume  fractions 
of  the  Si  active  sites,  b  was  kept  constant  at  1000  p,m,  while  a  was 
allowed  to  attain  values  of  lOOnm,  216  nm  and  480  nm  so  as  to 
vary  the  volume  fraction.  In  Fig.  3,  in  addition  to  the  energy  release 
rates  for  various  volume  fractions,  the  fracture  energy  of  the  matrix 
is  also  plotted.  Therefore,  based  on  Griffith’s  criterion,  the  crack 
radius  at  which  fracture  will  stop  is  the  point  at  which  the  energy 
release  rate,  G,  during  crack  growth,  exceeds  the  fracture  energy  of 
the  matrix.  Based  on  Fig.  3  it  can  be  seen  that  as  the  volume  fraction 
of  the  active  site  decreases,  cracking  will  stop  at  shorter  distances.2 


This  implies  that  anodes  containing  low  active  site  volume  fractions 
experience  less  cracking  and  therefore  the  anode  connectivity  is  not 
significantly  disrupted  during  cycling,  allowing  hence  for  the  mate¬ 
rial  to  respond  to  the  applied  voltages  during  cycling  and,  hence, 
providing,  therefore,  a  stable  capacity  retention. 

Despite  the  purely  mechanical  framework  employed,  the  afore¬ 
mentioned  predictions  by  Aifantis  et  al.  [20]  are  qualitatively  in 
agreement  with  the  experimental  evidence  shown  in  Fig.  1,  in 
which  it  is  clearly  seen  that  the  lowest  Si  content  anode  allowed  for 
a  more  stable  capacity  retention  than  the  higher  Si  content  anodes. 
More  mathematical  details  of  the  model  presented  above  can  be 
found  in  [18,20]. 


4.  Aspect  ratio  of  S11O2  particles  S11O2-C  nanocomposites 
anodes 


In  this  section  new  experimental  evidence  will  be  presented, 
that  documents  the  Sn  aspect  ratio  effect  in  capacity  retention.  This 
will  be  done  by  examining  the  microstructure  of  Sn/C  nanopowders 
that  can  be  used  as  anodes.3 

Fig.  4a  and  b  depicts  Sn  particles  attached  on  multiwalled  carbon 
nanotubes  [27].  It  can  be  seen  that  the  large  Sn  particles  covered 
the  whole  nanotube  surface,  forming  a  type  of  thin  film  layer.  Upon 
the  first  Li-insertion  the  capacity  was  above  500  mAh  g-1,  how¬ 
ever,  during  the  second  charge  the  capacity  was  negligible.  It  was 
anticipated  that  this  significant  capacity  drop  was  due  to  fracture 
and,  therefore,  scanning  (SEM)  and  transmission  (TEM)  electron 
microscopy  was  performed  after  the  first  cycle.  It  was  indeed  seen 
that  significant  cracks  formed  on  the  Sn  layer  (Fig.  4c),  while  fur¬ 
ther  TEM  indicated  that  the  large  Sn  particles  had  fractured  into 
multiple  nanoparticles  (Fig.  4d). 

In  addition  to  depositing  Sn  in  the  form  of  a  thin  layer  on  car¬ 
bon,  it  is  also  possible  to  attach  Sn  islands  on  amorphous  carbon 
obtaining  a  microstructure  of  fine  active  particles  dispersed  on 
the  carbon  surface.  For  this  case,  the  Sn  attaches  in  the  form  of 
Sn02  [27],  and  thus  far  a  fabrication  method  that  allows  the  attach¬ 
ment  of  pure  Sn  on  amorphous  carbon  does  not  exist.  Aifantis  et  al. 
[9]  cycled  the  material  of  Fig.  5  and  observed  a  stable  capacity  of 
385  mAh  g-1  for  500  cycles.  It  should  be  noted  that  the  capacity 
of  the  Vulcan-X72  that  was  used  as  the  matrix  was  180mAhg-1, 
prior  to  the  attachment  of  the  Sn02.  A  85%  capacity  increase,  hence, 
resulted  by  the  addition  of  8  wt%Sn.  Such  an  increase  in  capacity 
through  the  addition  of  such  a  low  wt%Sn  has  not  been  observed 
before;  note  for  example  that  in  [8]  50%  Sn  metal  was  added  to 
increase  the  graphite  capacity  by  40%.  The  electrochemical  stability 
observed  for  the  material  in  Fig.  5  indicates  that  the  Sn  nanoparti¬ 
cles  did  not  experience  fracture  upon  continuous  Li-insertion  and 
de-insertion,  similar  to  that  which  occured  for  the  Sn  film  layer  of 
Fig.  4. 

The  mechanics  based  advantage  of  the  island  design  (Fig.  5)  rel¬ 
ative  to  the  continuous  thin  film  geometry  (Fig.  4)  is  related  to 
the  resistance  of  the  Sn  island  particles  to  delamination  during 
volume  expansion.  This  resistance  to  delamination  of  islands  com¬ 
pared  to  thin  films  arises  from  the  stress  reduction  associated  with 
free  boundary  conditions  and  is  accentuated  at  a  high  aspect  ratio 
(island  height)/(island  width),  or  HI L.  For  a  given  strain  due  to  island 
volume  expansion  on  an  infinitely  stiff  substrate,  it  is  proposed  that 
the  critical  delamination  compressive  stress  concept  developed  in 
[28,29]  may  be  adapted  to  the  island  geometry  as 


(15)E  (2H\ 

c  12(l-i Mw 


(3) 


2  The  fracture  energy  of  C  could  not  be  found  in  the  literature  and  therefore  Y203 
was  used  as  the  matrix;  the  same  qualitative  results,  however,  would  be  true  for  - 

any  matrix  material.  3  These  materials  were  fabricated  and  presented  in  [27]. 
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Fig.  4.  (a)  SEM  image  of  as  prepared  Sn  attached  on  multi-walled  carbon  nanotubes;  Sn  appears  as  the  bright  phase,  (b)  TEM  image  of  as  prepared  Sn  attached  on  multi-walled 
carbon  nanotubes;  Sn  appears  as  black  phase,  (c)  SEM  image  after  first  de-lithiation  indicating  macro-cracks;  (d)  TEM  image  after  the  first  de-lithiation,  indicating  nanoscale 
fractured  particles.  SEM  scale  bars  are  10  p,m. 


where  H  is  the  island  half  thickness  and  L  is  the  island  half  length. 
Therefore,  if  the  compressive  stress  in  the  island  exceeds  the  critical 
delamination  stress,  then  the  island  will  detach  and  most  likely  fail 
to  contribute  to  the  electrochemical  activity.  It  can,  therefore,  be 


Fig.  5.  8  wt%Sn-Vulcan  C  nanopowder  prior  to  electrochemical  cycling.  Dark  areas 
indicate  Sn  rich  areas.  Representative  Sn  islands  are  contained  in  circles. 


understood  that  the  aspect  ratio  (H/L)  is  critical  to  the  mechanical 
stability  of  the  island  geometry.  In  order  to  test  the  idea  that  the  free 
surface  boundary  conditions  associated  with  the  island  geometry 
will  reduce  the  compressive  stress  relative  to  that  in  a  continu¬ 
ous  thin  film,  a  two  dimensional  elasticity  solution  is  pursued  for 
the  geometry  shown  in  Fig.  6.  The  stresses  and  displacements  that 
the  island  experiences  are  taken  to  be  given  by  the  biharmonic 
Marguerre  solutions  [30]  as 


crxx(x,y)= 


E 

T^V2 


f  d3V  d3V\ 
y  dx2dy+V  dy3  ) 


y(x,y)= 


1  — V2 


,  d3V  d3V 
t  +V  dx2dy+  dy3 


)■ 


<?xy(x,  y) 


E  f  d3#  d3#  \ 

1  -v2  y  dx3  +  dxdy2  J 


ux(x,y)=  — 


(4) 


where  the  energy  functional  ^  is  of  the  form 
^  =  cos(yftnX)  {A  sin  h{pny )  +  B  cos  h(/3ny)  +  C/3ny  sin  h(£„y));  (5) 
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Fig.  6.  Schematical  representation  of  Sn /C  island  microstructure. 
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Fig.  7.  Comparing  the  normalized  compressive  stress  for  different  H/L  values. 


where,  /3n  =  njrl2L  (n  =  odd),  while  the  constants  A ,  B,  are  found  from 
the  boundary  conditions 


—  0,  CTyy(X,  0)  —  0,  (Jxy(x,  0)  —  0, 


(6) 


as 


A  =  0  and  B  =  - 


2vC 
1  +  v 


C  is  found  from  the  definition  of  the  displacement  at  the 
islandicarbon  interface  as  a  linear  function  of  x  so  that 
A x/L  =  ux(y  =  -  2 H).  Expanding  the  displacement  as  a  Fourier  series 
in  which  C  becomes  a  Fourier  coefficient  (and  thus  a  function  of  the 
Fourier  index,  n,  the  various  values  of  Cn  may  be  determined  from 
Eq.  (7)): 


J  (^  sin (/8„x))  dx  =  J  [(ux(x,  -2H))sin(0„x)]  dx 


(7) 


where  A,  as  in  the  previous  sections,  denotes  the  unconstrained 
expansion  of  the  Sn  island. 

Hence, 


Cn  = 


n3n3 


32L3(v-  1)A  sin(n7r/2) 

(n7r)  [njrH(l  +p)cos  h{mzH / L)  -  L(v  -  l)sin  h(njrH/L)] 


Insertion  of  the  constants  of  integration,  A ,  B  and  C,  in  Eq.  (5)  and 
then  in  Eq.  (4)  allows  the  determination  of  the  compressive  stress 
axx\  A  is  taken  to  be  a  proportional  function  of  the  strain  associated 
with  the  maximum  constraint  to  island  expansion,  £,  A  =  sL.  The 
form  of  the  spatial  variation  in  the  normalized  compressive  stress 
along  x  at  y  =  (-H)  for  two  different  H/L  values  at  a  given  island 
volume  is  shown  in  Fig.  7.  It  may  be  seen  that  the  average  and 
maximum  amplitude  of  the  compressive  stress  is  reduced  in  value 
when  H/L  is  large  as  compared  to  when  it  is  small. 

Application  of  the  criteria  for  delamination  described  in  Eq.  (3) 
to  this  elasticity  solution  is  carried  out  by  calculating  the  average 
compressive  stress  in  an  island  and  then  comparing  this  aver¬ 
age  stress  to  the  critical  stress  for  delamination.  This  procedure 
is  developed  in  Fig.  8  for  an  island  structure  with  an  individual  Sn 
island  volume  of  16.75  nm3  is  obtained  by  calculating  the  average 
compressive  stress  in  an  island.  This  average  compressive  stress  is 
substituted  into  Eq.  (3)  giving 


2  L 

8HL2 


f  f 


axxdxdy  = 


15  E  (2  HY 


12(1  —  f2)  V  L 


(8) 


Fig.  8.  In  order  for  the  islands  to  remain  attached  to  the  carbon,  their  height/length 
ratio  (H/L)  must  be  greater  than  0.3. 

Rearrangement  of  terms  results  in  the  cancellation  of  elastic 
constants  or 


24(1  —v2)L 
32(15  )H3E  J 


f  f 


axxdxdy  =  1 


(9) 


Eq.  (9)  predicts  that  the  normalized  average  stress  must  be  less 
than  1  in  order  for  the  particle  to  be  attached  on  the  substrate  sur¬ 
face.  Plotting  therefore  the  left  hand  side  of  Eq.  (9)  as  a  function  of 
H/L,  it  is  found  that  it  intersects  the  constant  1  at  an  aspect  ratio  of 
approximately  0.3,  as  indicated  by  Fig.  8.  It  should  be  noted  that  the 
elastic  constants  of  Eq.  (9)  algebraically  cancel  and  do  not  influence 
the  solution. 

The  implication  here  is  that  for  a  given  percentage  of  volume 
change  during  cycling,  an  island  structure  has  a  greater  mechani¬ 
cal  stability  than  a  thin  film,  with  H/L  =  0.3  being  the  stability  limit; 
aspect  ratios  above  this  limit  increase  the  stability  of  island  attach¬ 
ment  on  the  substrate  during  cycling.  This  is  consistent  with  the 
evidence  of  Fig.  4,  which  shows  severe  fracture  of  a  Sn  thin  film 
during  the  first  cycle.  For  the  case  where  8  wt%  Sn  is  attached  in 
the  form  of  Sn02  on  Vulcan  carbon,  the  representative  Sn02  island 
in  the  white  circle  in  Fig.  5  has  a  height  of  2.86  nm  and  a  width  of 
5.2  nm,  which  indeed  gives  a  H/L  ratio  of  0.55,  which  is  the  above 
theoretical  limit  that  will  inhibit  fracture.  Furthermore,  the  islands 
in  black  circles  have  a  H/L  =  1 ,  which  again  inhibits  fracture.  There¬ 
fore,  fracture  was  avoided  for  this  material  and  a  high  capacity 
retention  was  obtained  [9]. 

5.  Conclusions 

In  the  present  study  it  was  illustrated  that  continuum  mechanics 
formulations  can  be  adapted  to  examine  fracture  and  damage  evo¬ 
lution  in  next-generation  anodes  for  rechargeable  Li-batteries.  This 
is  important  as  the  electrochemical  performance  of  these  power 
sources  is  directly  related  to  their  mechanical  stability.  In  the  first 
part  it  was  shown  that  according  to  Fig.  3  lower  Si  volume  frac¬ 
tions  embedded  in  a  less  active  matrix  experience  less  cracking  and 
therefore  the  anode  maintains  a  higher  mechanical  integrity  and  a 
better  capacity  retention  would  be  expected.  This  was  shown  to  be 
in  agreement  with  Si-SGG  nanocomposite  anodes  [7],  as  indeed  the 
anode  with  the  lowest  Si  volume  fraction  showed  better  capacity 
retention  than  the  anodes  containing  higher  Si  volume  fractions. 

In  the  second  part  of  this  paper  experimental  evidence  on  the 
mechanical  failure  of  Sn  thin  film  anodes  motivated  the  study  of 
the  mechanical  stability  in  nanoscale  island  structures.  Application 
of  linear  elasticity  indicates  that  the  aspect  ratio  of  the  Sn  islands 
plays  an  important  role  in  maintaining  the  mechanical  integrity  of 
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the  anode  during  cycling.  The  theoretical  prediction  reveals  that 
retaining  a  H/L  ratio  greater  than  0.3  will  inhibit  detachment  of  the 
Sn  islands.  This  was  found  to  be  in  agreement  with  the  microstruc¬ 
ture  of  a  SnO 2/C  anode,  whose  Sn02  islands  had  a  H/L  ~  0.5-1  and 
gave  a  stable  capacity  for  over  300  cycles  [9],  whereas  a  contin¬ 
uous  Sn  layer  fractured  severely  after  the  first  cycle  as  shown  in 
Fig.  4.  It  should  be  noted  that  other  studies  (e.g.  [31])  have  also 
indicated  that  Sn  based  thin  films  exhibit  a  poor  capacity  retention 
and  attributed  it  to  fracture.  However,  TEM  and  SEM  images  docu¬ 
menting  this  during  the  first  cycle  or  a  detailed  theoretical  analysis 
describing  the  underlying  damage  mechanics  have  not  been  pre¬ 
sented  before. 

In  concluding,  it  should  be  noted  that  in  addition  to  volume  size 
and  aspect  ratio  of  the  particles,  particle  size  also  plays  a  significant 
role  in  fracture,  as  first  noted  by  [32].  Such  theoretical  considera¬ 
tions  are  currently  being  undertaken. 
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